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The characteristics of Fe-based superconductors are manifested in their electronic, magnetic prop¬ 
erties, and pairing symmetry of the Cooper pair, but the latter remain to be explored. Usually in 
these materials, superconductivity coexists and competes with magnetic order, giving unconven¬ 
tional pairing mechanisms. We report on the results of the bulk magnetization measurements in 
the superconducting state and the low-temperature specific heat down to 0.4 K for BaFe 2 -a:Nij; As 2 
single crystals. The electronic specihc heat displays a pronounced anomaly at the superconduct¬ 
ing transition temperature and a small residual part at low temperatures in the superconducting 
state. The normal-state Sommerfeld coefficient increases with Ni doping for x = 0.092, 0.096, and 
0.10, which illustrates the competition between magnetism and superconductivity. Our analysis of 
the temperature dependence of the superconducting-state specihc heat and the London penetration 
depth provides strong evidence for a two-band s-wave order parameter. Further, the data of the 
London penetration depth calculated from the lower critical held follow an exponential temperature 
dependence, characteristic of a fully gapped superconductor. These observations clearly show that 
the superconducting gap in the nearly optimally doped compounds is nodeless. 

PACS numbers: 74.25.Bt, 74.25.Dw, 74.25.Jb, 65.40.Ba 


I. INTRODUCTION 


One of the major themes in the physics of condensed 
matter is unconventional superconductivity in Fe-based 
materials^^'*’. These materials have multiple Fermi pock¬ 
ets with electronlike and holelike dispersion of carriers 
and both hole and electron Fermi pockets show a low car¬ 
rier density®. Superconductivity appears at the border of 
the antiferromagnetic (AF) regime, which may have a sig¬ 
nificant impact on the pairing mechanism®. However, the 
exact picture of the interplay between superconductivity 
and magnetism remains elusive^. Although other sce¬ 
narios involving orbital fluctuations are possible, it has 
generally been believed that spin fluctuations play an im¬ 
portant role and act as the mediating bosons for electron 
pairing and superconductivity'. Despite great successes 
in studying these materials, there are still unresolved is¬ 
sues, particularly the symmetry and structure of the or¬ 
der parameter, and doping evolution of the supercon¬ 
ducting (SC) gap, which should provide an understand¬ 
ing of the pairing mechanism of these systemsIt has 
been well characterized that both cuprates and conven¬ 
tional phonon-mediated superconductors are character¬ 
ized by distinct d-wave and s-wave pairing symmetries 
with nodal and nodeless gap distributions, respectively. 


There is no general consensus on the nature of pairing 
in iron-based superconductors leaving the perspectives 
ranging from 5 '“'"''" wave, to S^, and to d wave®“^^. In 
addition, from ®®Co and As nuclear magnetic resonance 
measurements, the spin triplet order parameter was ruled 
out in BaFei.8Coo,2As2^®. It turns out that the SC gap 
distributions are vary with different systems and unusu¬ 
ally are sensitive to the sample quality. For systems 
with both hole and electron Fermi surfaces, such as op¬ 
timally doped Bai_2,K2;Fe2As2®^, Ba(Fei_a;Coa;)2As2^®, 
NaFei_2,Coa;As^®, and Fe(Se, Te)^®, the gaps mea¬ 
sured by low temperatures specific-heat fit well to the 
predictions of two nodeless SC gap. The tempera¬ 
ture dependence of the lower critical field in LiFeAs^®, 
Bao. 6 Ko. 4 Fe 2 As 2 ^^, FeSe^®, and Cao. 32 Nao. 68 Fe 2 As 2 ® has 
supported the existence of two s-wave-like gaps. The pos¬ 
sibility of nodes along the c axis in the superconducting 
gap has been reported in NdFeAsO0.82F0.i8 and LaFePO, 
where the magnetic penetration depth exhibited a nearly 
linear temperature dependence^®’®®. 

We begin with listing several facts about 
BaFe2_xNia;As2. (i) In the first, it appears as an 
ideal candidate to study the fundamental properties of 
superconductivity due to the availability of high-quality 
single crystals with rather large dimensions®®, (ii) In the 
undoped state, BaFe2As2 shows a combined spin-density 
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wave (SDW) and structural transition near Tjv = Tg 
= 138 K. The pristine compound is characterized by a 
bad metallic behavior with a coherent Drude component 
and doping with Ni and P transforms a bad metal to 
a good metal, while the system remains a bad metal 
with K-doping^^. (iii) The Neel temperature of the 
electron-doped iron pnictides decreases gradually with 
increasing electron-doping level, and the AF phase 
appears to coexist with the SC phase^^“^®. However, a 
neutron-scattering study reveals an avoided quantum 
critical point, which is expected to influence the proper¬ 
ties of both the normal and SC states strongly^^. This 
raises a critical question concerning the role of quantum 
criticality^® and the coexistence of magnetism and 
superconductivity to the SC pairing structure®^, (iv) 
Furthermore, a recent neutron-scattering measurement 
has revealed that the low-energy spin excitations in 
BaFe 2 -a;Ni 2 ;As 2 change from fourfold symmetric to 
twofold symmetric at temperatures corresponding to 
the onset of the in-plane resistivity anisotropy. In the 
overdoped compounds both resistivity and spin excita¬ 
tion anisotropies are vanished. Therefore, they are likely 
intimately connected®, (vi) The London penetration 
depth A measurements suggest that the competition 
between superconductivity and magnetic/nematic order 
in hole-doped compounds is weaker than in electron- 
doped compounds®^. In this context, it is important to 
understand the doping, field, and temperature depen¬ 
dence of AF spin correlations. Studying the symmetry 
and structure of the order parameter is a key not only 
to understand all these interesting features but also to 
address unsettled issues in BaFe 2 -a;Nia;As 2 . 

Low-temperature specific heat Cp and the London 
penetration depth A are two powerful techniques for 
probing the gap structure of bulk superconductors. Both 
measurements probe bulk SC properties. A is a funda¬ 
mental parameter which detects the pairing symmetry 
and the T-dependence of A can determine gap function. 
Since Cp is directly related to the quasiparticle density 
of states, its temperature dependence reflects the nature 
of the SC state such as gap symmetry, the presence of 
multigaps, and coupling strength between electrons and 
phonons. In addition, it is less affected by vortex pin¬ 
ning. An exponential vanishing of the specific heat at 
low temperature in conventional s-wave superconductors 
is caused by the finite gap in the quasiparticle spectrum. 
This is due to the quasiparticle thermal fluctuations go 
exponentially to zero as T —>• 0. For SC gap with gap 
nodes, electronic excitations are possible even at very 
low temperatures^^. In general, specific heat comprises 
of two parts: the electronic Cei and the phononic contri¬ 
bution Cph- Information about the pairing symmetry is 
contained in the C^i , which is proportional to the density 
of states (DOS) at the Fermi energy. Exploring the sym¬ 
metry and structure of the order parameter, and the evo¬ 
lution of the SC gap with Ni doping in BaFe 2 -a;Nia;As 2 
system based on the mentioned two bulk detection tech¬ 
niques is thus highly desired. It should be mentioning 


that we have estimated the Cph from BaFei, 75 Nio. 25 As 2 . 
This sample is not superconductor throughout the tem¬ 
perature range as evident in Fig. 2(e) where the C/T ex¬ 
hibits a monotonous increase against the temperature. 
The fact that the low temperature specific heat data for 
the investigated samples [the inset of Fig. 2(e)] exhibit 
a linear behavior at low temperatures without any up¬ 
turn indicates the absence of Schottky-like contributions 
in our samples. Furthermore, at T > Tc the specihc heat 
data of the SC and non-SC samples are comparable, con¬ 
firming similar phonon contributions to the specific heat 
of SC and non-SC samples. Therefore, magnetic contri¬ 
bution to specific heat will be negligible and the specihc 
heat can be assumed to have contribution from the elec¬ 
tronic and lattice part only. 

In this work, we use magnetization and low temper¬ 
ature specihc heat measurements on BaFe 2 -a;Nia;As 2 to 
study the interplay between magnetism and supercon¬ 
ductivity with the emphasis of nature of the SC pairing 
symmetry by focusing on materials near optimal doping 
[Fig. 1(a)]. Based on the comprehensive low-T measure¬ 
ments, we provide evidence for nodeless superconductiv¬ 
ity in the doping range of x = 0.092, 0.096, and 0.10. 
The temperature-dependence of Ao&(T) calculated from 
the lower critical held and the Cei can be well described 
by using a two-band model with s-wave-like gaps. Re¬ 
liable values of the normal-state Sommerfeld coefficients 
are obtained for the studied materials, which increases 
with Ni doping, illustrating the strong competition be¬ 
tween magnetism and superconductivity. 


II. EXPERIMENTAL 

BaFe2-:rNRAs2 (x = 0, 0.03, 0.065, 0.092, 0.096, 0.10, 
0.15, and 0.25) single crystals were grown by the FeAs 
self-hux method®^. The actual Ni level was determined 
to be 80% of the nominal level x through the inductively 
coupled plasma analysis of the as-grown single crystals. 
Magnetization measurements were performed by using a 
Quantum Design SC quantum interference magnetome¬ 
ter. The low-T specihc heat down to 0.4 K was measured 
in its Physical Property Measurement System with the 
adiabatic thermal relaxation technique along H \\ c up to 
H = 9T. 


III. RESULTS AND DISCUSSION 

The arrows in Fig. 1(a) indicate eight doping levels 
presented in this work. These include x = 0 (parent 
compound shows a Ts{Tm) = 137(2) K), x = 0.03 and 
0.065 (lightly electron-doped non SC and SC samples 
with Ts/Tm = 110/104 and 82/70 K) respectively, x = 
0.092, 0.096, and 0.10 (nearly optimal doping SC samples 
with static incommensurate short-range order), x = 0.15 
(overdoped superconducting sample without AF order 
coexisting with superconductivity) and x = 0.25 (heav- 
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FIG. 1: (a) The electronic phase diagram of BaFe 2 - 2 ;Nia;As 2 obtained from magnetic and specific heat data, showing the 

suppression of the magnetic (Tn) and structural (Tg) phase transitions with increasing Ni concentration and the appearance of 
the SC transitions. The arrows indicate eight doping levels studied in this work. The PM Tet, PM Orb, C-AF, and IC-AF are 
paramagnetic tetragonal, paramagnetic orthorhombic, commensurate AF orthorhombic, and incommensurate AF orthorhombic 
phases, respectively. The inset illustrates the electron-doping dependence of Tg — Tn- (b) shows the temperature dependence of 
the magnetic susceptibility in an external field of 10 Oe applied along the c axis. The susceptibility has been deduced from the 
dc magnetization measured by following ZFC and FC protocols of BaFe 2 -a:Ni 2 ,As 2 single crystals, (c) presents the isothermal 
magnetization M vs. H loops measured at 2 K up to 9 T for H \\ c for x = 0.092, 0.096, 0.10, and 0.15. (d) 14-18.5 K for each 
0.5 K, (e) 13-18.5 K for each 0.5 K, (f) 13-19 K for each 0.5 K, and (g) 7-13 K for each 0.5 K plots at high temperatures exhibit 
a pronounced second peak for x = 0.092, 0.096, 0.10, and 0.15 respectively. 


ily overdoped non SC sample). As for BaFe 2 -a;Coa;As 2 
system^®, it has been shown that near optimal super¬ 
conductivity (see Fig. 1(a)), the commensurate static AF 
order changes into transversely incommensurate short- 
range AF order that coexists and competes with super¬ 
conductivity®®. Similar to the case of BaFe 2 _a;Coa;As 2 
and CaFe 2 -a;Coa;As 2 , the underdoped region exhibits a 
splitting of the structural and magnetic phase transitions. 
The inset of Fig. 1(a) shows the electron-doping depen¬ 
dence of Ts—Tn- Figure 1(b) shows the magnetic suscep¬ 
tibility measured with the zero field cooling (ZFC) and 
field cooling (FC) in an external field of 10 Oe applied 
along the c axis. The FC and ZFC data prove a sharp dia¬ 
magnetic signal. Beyond, the SC volume fraction is close 
to 1, thus confirming bulk superconductivity and the high 
quality of BaFe 2 -xNia:As 2 single crystals. The Tc has 
been determined from the onset diamagnetic transition 
temperature between ZFC and FC to be around ^ 7.7, 
18.5, 19.0, 20.0, and 13.9 K for x = 0.065, 0.092, 0.096, 
0.10, and 0.15 respectively. The clear irreversibility be¬ 
tween FC and ZFC measurements is the consequence of 
a strong vortex trapping mechanism, either by surface 
barriers or bulk pinning. 

Figure 1(c) presents the field dependence of the isother¬ 
mal magnetization M at 2 K up to 9 T for i7 || c for x 
= 0.092, 0.096, 0.10, and 0.15. At T = 2K for x = 
0.096 and 0.10, the M{H) exhibits irregular jumps close 
to H = 0 similarly to LiFeAs, Bao. 65 Nao, 35 Fe 2 As 2 , and 
Cao, 32 Nao. 68 Fe 2 As 2 superconductors®’^®’^®. Figures l(d- 
g) present the field dependence of the isothermal mag¬ 
netization M at various temperatures very close to 


up to 9T for X = 0.092, 0.096, 0.10 and 0.15 respec¬ 
tively. In addition, the SC M{H) exhibits no magnetic 
background. This indicates that our investigated sam¬ 
ples contain negligible magnetic impurities. The width 
of the magnetic loops decrease while increasing the ap¬ 
plied field. However, at higher temperatures the width of 
the loops initially decreases showing a minimum at the 
Hm field and then increases again. Further, the M{H) 
loops demonstrate another pronounced peak or so-called 
second peak. The second peak effect has been studied 
extensively and its origin may be attributed to various 
mechanisms. It has been well established that the sec¬ 
ond peak effect is strongly influenced by the oxygen de¬ 
ficiency in cuprates'®®’^®. In the case of Fe-based super¬ 
conductors, the local magnetic moments may form the 
small size normal cores, and may be a possible reason 
of the second peak effectHowever, the real pinning 
mechanism needs further investigation. The position of 
the second peak shifts to higher fields while decreasing 
temperature, eventually beyond the available field range. 
This can explain the nonvisibility of a second peak at low 
temperatures in the Figs, l(d-g). The M{H) loops show 
irreversibility in magnetization, which vanishes above a 
characteristic field Hir-r [Figs, l(d-g)]. It is noteworthy 
that the first vortex penetration field may not reflect the 
true Hci{T) because of Bean-Livingston surface barrier. 
The fact that the hysteresis loops for x = 0.092, 0.096, 
O.IO, and 0.15 are symmetric around M = 0, pointing 
to relatively no surface barriers and implying that the 
bulk pinning plays a dominant role in our investigated 
compounds. In contrast to that, if surface barriers were 
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FIG. 2: (a) Temperature dependence of the specific heat of BaFe 2 -xNi 3 ;As 2 (x = 0, 0.03, 0.065, 0.092, 0.096, 0.10, 0.15, and 

0.25) measured in zero magnetic field. The insets show the derivative of specific heat for the crystal with x = 0.03 (upper inset) 
and X = 0.065 (lower inset), where structural and SDW transitions can be clearly recognized from the dips, (b) x = 0.10, (c) 
X = 0.096, and (d) x = 0.092 indicated the neutron data counting time (30 min=point on HB-IA taken from Ref. [35]) and 
the specific heat. The Tjv is marked by an arrow, (e) The temperature dependence of the specific heat C/T of samples with x 
= 0.092, 0.096, 0.10, and 0.25 down to T = 400 mK. The upper inset shows the low-temperature specific heat of two samples 
with X = 0.10 and 0.25. The straight lines represent linear fits to Cp = ')T + The lower inset presents the enlarged CpfT 
vs. T plot near the SC transition for x = 0.096. The lines show how CpjT^ and T, are estimated. 


predominant, the first vortex entrance can occur at much 
higher field (« He). This is a very important point in or¬ 
der to obtain reliable estimations of the thermodynamic 
lower critical field (see below). 

Specific heat provides a probe for the symmetry and 
structure of the SC order parameter. Figure 2(a) sum¬ 
marizes the temperature dependence of the zero-field 
specific-heat data at various Ni-doping levels in the 
BaFe 2 -a:NixAs 2 series plotted as Cp vs. T. The data 
of the parent-compound {x = 0) shows a very sharp 
first-order structural transition coinciding with the SDW 
transition at 136 K (upon heating) and with a transi¬ 
tion width of about 3K. Because of the narrowness of 
the transition, a temperature rise of only 0.5 % was used 
for each measurement in the vicinity of the transition 
of all measurements. Upon Ni-doping, the sharp first- 
order structural/magnetic anomaly of the parent com¬ 
pound gradually broadens, shifts and splits to lower tem¬ 
peratures and is considerably reduced in magnitude. For 
X = 0.03 and 0.065, the combined structural/magnetic 
anomaly of the pristine compound actually splits into 
two distinct anomalies at 110, 104 and 86, 74 K, respec¬ 
tively. The error in the determination of the Tg and Tjv 
transition temperatures can be estimated at around 2 K 
if we take into account that the peak in the first deriva¬ 
tive of the specific heat is relatively sharp [see upper and 
lower insets of Fig. 2(a)]. Then, the transition is shifted 
to 40(4), 30(5), and 32(5) K for x = 0.092, 0.096, and 
0.10, respectively [see Fig. 2(b-d)]. These data are in line 
with the recent high-resolution x-ray and neutron scat¬ 
tering data as discussed in Ref. [35]. Recent neutron 
scattering data on x = 0.10 sample reveal a weak static 
AF order with magnetic scattering 5 times smaller than 
that of X = 0.096. In spite of the small moments of x 
= 0.10, the temperature dependence of the magnetic or¬ 


der parameters for both samples indicates that their AF 
temperatures are essentially unchanged at ± 5K [35]. 

Figure 2(e) shows the temperature dependence of the 
specific heat of the samples with x = 0.092, 0.096, 0.10, 
and 0.25 down to 0.4 mK. An entropy conserving con¬ 
struction has been used to determine the SC transition 
temperature from the specific heat data. For x = 0.092, 
0.096, and 0.10 a clear anomaly at 18.4, 18.9, 20 K respec¬ 
tively indicates the onset of bulk superconductivity. The 
sample, with x = 0.25, remains in the normal state. The 
fact that the low temperature specific heat data exhibit 
a linear behavior at low temperatures without any up¬ 
turn indicates the absence of Schottky-like contributions 
in our investigated samples [see upper inset of Fig. 2(e)]. 
It is important to note that it is impossible to obtain 
the lattice background by fitting the specific heat of the 
SC samples to an odd-power polynomial above Tc due to 
the electronic term of the total signal of the specific heat 
data. As demonstrated below, a more reliable phonon 
term can be estimated from the data of the x = 0.25 
sample, whose low-temperature specific heat follows pre¬ 
cisely the Debye law between 0.4 and 4.5 K, with 7 „ = 
22.9mJ/mol and (3 = 0.23mJ/mol K'*. 

Further experimental investigations on the structure 
and magnitude of the SC gaps in BaFe 2 -xNixAs 2 by 
means of bulk specific heat data are of great interest. 
In order to determine the specific heat related to the 
SC phase transition we need to estimate the Cph and 
Cel contributions to Cp in the normal state. In or¬ 
der to determine the phononic contribution to the spe¬ 
cific heat for x = 0.25, the following relation is used: 
(7 x=0.25 ^ c^=o-^5 _ Cff is just yT. 

The same shape of the phononic heat capacity in the SC 
samples and overdoped sample is assumed. Therefore, 
the specific heat of the SC samples can be represented 
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by: 

= Cfi/T - (1) 

which allows us to calculate the Ce\ of the SC samples. 
The small deviation of the scaling factor g from unity, 
plausibly related to experimental uncertainties, demon¬ 
strates that the above procedure represents a very good 
method to determine the phonon background. The value 
of g was determined from the requirement of equality be¬ 
tween the normal and SC state entropies at Tc, that is 
(Cei/T) dT = ynTc, where 7 n is the normal state elec¬ 
tronic specific heat coefficient. We started with g = 1, 
but we found that the entropy conservation criterion is 
satished with g = 0.95. Physically, this indicates that 
the substitution of Fe by Ni does not substantially affect 
the lattice properties. 

Figure 3 shows the temperature dependence of the elec¬ 
tronic contribution to the specific heat in the zero field 
determined by subtracting Cph for x = 0.10 [Fig. 3(a)], 
0.096 [Fig. 3(b)] and 0.092 [Fig. 3(c)]. The entropy con¬ 
servation required for a second-order phase transition is 
fulfilled as shown in the inset of Fig. 3(b). This check 
warrants the thermodynamic consistency for both, the 
measured data and the determination of Cei. It is obvi¬ 
ous from Fig. 3 that the SC transition at Tc is well pro¬ 
nounced showing a sharp jump in Cgi at Tc. The jump 
height of the specific heat at Tc is found to be ACei jT^ ~ 
23(0.5), 24.8(2), and 25.1(1) mJ/mol for a; = 0.092, 
0.096, and 0.10, respectively. Generally, the specific heat 
jumps at Tc obtained for these materials scale relatively 
well with its Tc in light of the recent careful results for 
the pnictide superconductors^’^^ in which the universal 
curve ACp/Tc oc T^ is explained. Furthermore, it has 
been well reported that the jump of the specific heat 
AC/Tc varies with Tc, and has a peak near optimal dop¬ 
ing and decreases at smaller and larger doping. This is 
a direct manifestation of the coexistence between anti¬ 
ferromagnetism and SC order parameters'^®. From our 
determined 7 n = 20.5, 23.5, and 24.6mJ/mol for x 
= 0.092, 0.096, and 0.10 respectively, we find ACei/ynTc 
= 1.1, 1.06, and 1.04 for x = 0.092, 0.096, and 0.10 re¬ 
spectively. These values are smaller than the prediction 
of the weak coupling BCS theory (ACei/ 7 nTc = 1.43). 
Taking into account the fact that the SC transition is 
relatively sharp in our SC samples, a distribution in Tc 
or the presence of impurity phases cannot explain the 
reduced value of the specific heat jump. In addition, 7 n 
increases with Ni doping, illustrating the competition be¬ 
tween magnetism and superconductivity. 

We believe, however, that the presence of multiple SC 
gaps may reduce the universal parameter, as evidenced 
in other 122 Fe-based superconductors^^. It has been 
also well reported that the reduced jump in the specific 
heat ACp/Tc compared to that of a single-band s-wave 
superconductor might be related to a pronounced multi¬ 
band character with rather different partial densities of 
states and gap values^"^. Note that Cei/T almost satu¬ 
rates at low temperature; however, it does not extrap- 



FIG. 3: The electronic specific heat Cei/T as a function of 
temperature for BaFe 2 -a;Nia;As 2 [a: = 0.1 (a), 0.096 (b), and 
0.092 (c)]. The inset in (b) presents the entropy in the normal 
and superconducting state as a function of T. The inset in 
(c) shows the low-T data on a larger scale. 7n represents 
the normal-state electronic coefficient of the specific heat and 
7 r is the residual electronic specific heat. The dashed lines 
represent the theoretical curves based on single-band weak- 
coupling BCS theory, while the solid lines illustrate the d- 
wave approximation. The solid red lines indicate the curves 
of the two s-wave gap model. 


olate to zero, yielding a residual electronic specific-heat 
value 7 r = 2.6, 0.9, and I.6mJ/mol for x = 0.092, 
0.096, and O.IO, respectively. The finite value of 7 r indi¬ 
cates a finite electronic density of states at low energy, 
even in zero applied field. We mention that the pres¬ 
ence of a finite 7 r is common in both electron- and hole- 
doped 122 crystals and that the value of 7 r is remarkably 
low, showing the good quality of our investigated sin¬ 
gle crystals. However, the origin of this residual term 
is still unclear. It may be because of an incomplete 
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transition to the SC state or because of broken pairs 
caused by disorder or impurities in unconventional su¬ 
perconductors, and/or spin-glass behavior. On the other 
hand, previous specific heat measurements on optimally 
doped YBa 2 Cu 3 07_5 exhibit such a 7 r term. For in¬ 
stance, even the best YBa 2 Cu 3 06.56 samples present 7 r 
« 1.85 mJ/mol [48]. It has been proposed that this 7 ^ 
term originates from a disorder-generated finite density 
of quasiparticle states near the d-wave nodes. It is worth 
to mention that 7 r in our SC samples reaches 12.6, 3.8, 
and 6.5% of the normal state Sommerfeld coefficient 7 n 
for X = 0.092, 0.096, and 0.10, respectively. A similar ob¬ 
servation of (7r/7n ~ 5.7%- 24 %)^"^ d9 ^as also reported 
in iron pnictide superconductors. 

The almost linear temperature dependence of Cei/T 
of the SC samples indicates that the specific heat data 
cannot be described by a single BCS gap. In order to 
illustrate this we show a theoretical BCS curve with 
A = 1.764 ^bTc = 2.23 meV in Fig. 3. One can see that 
systematic deviations from the data are observed in the 
whole temperature range below T^. Since a single gap 
cannot describe the data, we applied a d-wave calculation 
and a phenomenological two-gap model developed for the 
specific heat of MgB 2 ®° as in Eq. (2) and Eq. (3). Eor the 
d-wave approximation we used A = Aocos(2d). In the 
case of a two-band model, the thermodynamic properties 
are obtained as the sum of the contributions from the in¬ 
dividual bands, i.e., ai = Ai/fc^Tc and a 2 = A 2 /kBTc 

c c A roo 

[/ln/ + (l-/)ln(l-/)]d,, (2) 
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(3) 


where t = T/Tc, f = [exp(/3£'-|- 1)]“^, /3 = and 

the energy of the quasiparticles is given by if = [e^ -|- 
A2(^)]0.5 g being the energy of the normal electrons 
relative to the Fermi level. The integration variable is 
y — c/Aq. In Eq. (2), the scaled gap a = Ao/ksT is 
the only adjustable fitting parameter in the case of a 
single gap. At the same time, 7i/7n (* = 1, 2), which 
measure the fraction of the total normal electron density 
of states, are introduced as adjustable parameters. This 
fitting is calculated as the sum of the contributions from 
two bands by assuming independent BCS temperature 
dependencies of the two SC gaps. 

The best description of the experimental data for each 
type of order parameter, d-wave and two-gaps s-wave can 
be seen in Fig. 3. More obvious deviations exist in the 
case of the d-wave approach for the SC samples. This 
clearly indicates that the gap structure of our systems 
is more likely to be nodeless s-wave, which is reasonably 
well comparable with the penetration depth data (see be¬ 
low). The good description of the experimental data for 
the two-gaps s-wave model is obtained by using values of 
Ai(0) = 1.74, 1.8, and 1.85 kgT^, A2(0) = 0.68, 0.74, and 
0.79 ksTc for x = 0.092, 0.096, and 0.10, respectively. For 


FIG. 4: The upper panel shows the phase diagram of Hci 
vs. the applied temperatures of BaFe 2 -a;Nia;As 2 (x = 0.092, 
0.096, and 0.10) for the field applied parallel to the c axis. 
The bars show the uncertainty of estimation by the deviat¬ 
ing point of the regression fits. The error bar in the values of 
Hci is about 5 Oe of the investigated samples. The upper inset 
shows the field dependence of the superconducting initial part 
of the magnetization curves measured of BaFei. 9 oNio.ioAs 2 at 
various temperatures for H || c. The lower inset depicts an 
example used to determine the Hd value using the regres¬ 
sion factor, R, at T = 2 K. The lower panels present the field 
dependence of the typical plot of \/Mt vs H at various tem¬ 
peratures for X = 0.10. The solid lines are a linear fit to the 
high-field data of y/Mt vs. H. Hd values are determined by 
extrapolating the linear fit to y/Mt = 0. 


the investigated systems, the large gap A^ has a higher 
value than the weak-coupling BCS (1.76fcBTc) gap value, 
while the smaller one As has a value lower than the BCS 
one. This is consistent with the theoretical constraints 
that one gap must be larger than the BCS gap and one 
smaller in a weakly coupled two-band superconductor®^. 
Similar studies have been outlined in iron-based super¬ 
conductors (Table I). 

Next we discuss the temperature dependence of the 
lower critical field Hd, the field at which vortices pene- 
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trate into the sample, in the SC-state, which is another 
independent test sensitive to the gap structure. However, 
determining the Hd from magnetization measurements 
has never been an easy task. In order to determine the 
exact values of the Hd from the low-field M-H curves 
measured at different temperatures, we have to detect the 
onset of the small deviation from the perfect diamagnetic 
signal. This is rather difficult and sometimes a debat¬ 
able process. The most popular method to estimate Hd 
consists of detecting the transition from a Meissner-like 
linear M{H) regime to a non-linear M{H) response (see 
the upper inset of Fig. 4 upper panel), once the vortices 
penetrate into the sample and build up a critical state. 
This transition is not abrupt therefore bearing a substan¬ 
tial error bar. These sort of measurements are obtained 
by tracking the virgin M{H) curve at low fields at sev¬ 
eral temperatures, as shown in Fig. 4 for H \\ c for x = 
0.092, 0.096, and 0.10. These M-H curves show at low H 
a linear dependence of magnetization on the field indica¬ 
tive of Meissner phase as well deviation from linearity 
at higher fields. We have adopted a rigorous procedure 
{i.e. user-independent outcome) to determine the transi¬ 
tion from linear to non-linear M{H), which consists of 
calculating the regression coefficient i? of a linear fit to 
the data points collected between 0 and H, as a func¬ 
tion of H. Then, Hd is taken as the point where the 
function R{H) departs from 1. The result of these cal¬ 
culations is illustrated in the lower inset of Fig. 4 upper 
panel. Additionally, the temperature dependence of the 
first vortex penetration field has been experimentally ob¬ 
tained by measured the onset of the trapped flux moment 
Mt as described in Refs. [52,53]. In contrast to tracking 
the virgin M{H) curves at low fields at several tempera¬ 
tures where a heavy data post-processing is needed now 
a careful measurement protocol needs to be followed with 
little data analysis. Indeed, the Hd values obtained from 
the onset of the Mt are close to those obtained from the 
latter method. 

Once the values of Hd have been experimentally de¬ 
termined, we need to correct them accounting for the de¬ 
magnetization effects. Indeed, the deflection of field lines 
around the sample leads to a more pronounced Meissner 
slope given by M/Ha = —^/{^ — N)j where N is the de¬ 
magnetization factor. Taking into account these effects, 
the absolute value of Hd can be estimated by using the 
relation proposed by Brandt®"^: 

qdisk = ^ + ^tanh[1.27- ln(l -h ^)], (4) 

6n Sty a o 

where q = {\M/Ha\ — l){b/a), and a is the average of 
the dimensions perpendicular to the field of our investi¬ 
gated sample. For our samples we find N « 0.958(0.1), 
0.95(0.12) and 0.94(0.085) for x = 0.092, 0.096, and 0.10, 
respectively. The corrected values of Hd obtained by fol¬ 
lowing the two methods described above, are illustrated 
in the main panel of Fig. 4 for H \\ c. In fact, the determi¬ 
nation of Hd allows us to extract the magnetic penetra¬ 
tion depth, a fundamental parameter characterizing the 



_I_ I _I_ I _I_ I _ 1 _ I \ \ _ 

0 4 8 12 16 20 

T(K) 


FIG. 5: The T-dependence of the Xab{T) for BaFe 2 _a:Nij,As 2 
[x = 0.10 (a), 0.096 (b), and 0.092 (c)]. The red solid lines 
are the fitting curves using a two-gap model. The solid and 
dashed lines represent the d-wave and a single-gap BCS ap¬ 
proach, respectively. The inset of (b) presents the tempera¬ 
ture dependence of the magnetic penetration depths Xab- 

SC condensate which carries information about the un¬ 
derlying pairing mechanism. In the SC state, the temper¬ 
ature dependence of the penetration depth is a sensitive 
measure of low-energy quasiparticles, making it to a pow¬ 
erful tool for probing the SC gap^®. In order to shed light 
on the pairing symmetry in our system, we estimated the 
penetration depth at low temperatures using the tradi¬ 
tional Ginzburg-Landau theory, where Hd is given by: 
^J-oHa/ ~ ('(’o/47rA^^) In Kc, where 4>q is the magnetic-flux 
quantum = h/e* = 2.07xl0”^0e cm^, Kc =^ab/^ab is 
the Ginzburg-Landau parameter^®’®®, which we obtained 
at A(0) =214(15), 255(10), and 240(10) nm for x = 0.092, 
0.096, and 0.10, respectively. 

The temperature dependence of the Xab applied along 
the c axis is shown in the inset of Fig. 5(b). At low 
temperatures from the inset, Xab(T) does not show an 
exponential behavior as one would expect for a fully 
gapped clean s-wave superconductor. The main fea¬ 
tures in Fig. 5, A(r)-data, can be described in the fol¬ 
lowing ways: (i) As the first step we compare our data 
to the d-wave and single-gap BCS theory under the 
weak-coupling approach (see solid and dashed lines in 
Fig. 5). Indeed, both quantities lead to a rather dif- 











TABLE I: The superconducting transition temperature Tc (in K), the SDW transition temperature Tjv (in K), the residual 
and normal-state electronic specific heat 7 r and 7 „, respectively (in mj/mol K^), the universal parameter ACei/7nTc, and the 
superconducting gap properties extracted from specific-heat and lower critical field (-ffci) measurements for BaFe 2 -i:Ni 2 ;As 2 
(x=0.10, 0.096, and 0.092) along with other 122 Fe-based superconductors. 


Compounds 

Tc 

Tn 

7r 

7n 

ACel/ynTc 

<1 

As/ksTc 

Al/As 


Technique 

Ref. 

BaFei.9oNio.ioAs2 

20(1) 

30(3) 

1.6 

24.6 

1.04 

1.85, 1.9 

0.79, 0.68 

2.3, 2.7 

0.41, 0.59 

C{T),Xcb 

this work 

BaFei.904Nio.o96As2 

19(0.5) 

32(5) 

0.9 

23.5 

1.06 

1.8, 1.74 

0.74, 0.59 

2.4, 2.9 

0.44, 0.56 

C{T),Xab 

this work 

BaFei,9ogNio.o92As2 

18.4(0.2) 39(4) 

2.6 

20.5 

1.12 

1.74, 1.72 

0.68, 0.49 

2.5, 2.9 

0.39, 0.61 

C{T),Xab 

this work 

Ba(Feo.925Coo.075)2As2 

21.4 

- 

5.77 

23.8 

1.2 

2.2 

0.95 

2.3 

0.33, 0.67 

C{T) 

[47] 

Bao.6Ko.4Fe2As2 

35.8 

- 

1.2 

50 

1.54 

2.88(0.2) 

0.64(0.02) 4.45(0.3) 

0.5, 0.5 

C{T),Hci 

[11,27] 

Bao.65Nao.35Fe2As2 

29.4 

- 

3.3 

57.5 

1.26 

2.08 

1.06 

1.96 

0.48, 0.52 

C{T) 

[49] 


ferent trend and show a systematic deviation from the 
data in the whole T-range below T^. (ii) Then, the ob¬ 
tained temperature dependence of was analyzed 

by using the phenomenological a-model. This model 
generalizes the temperature dependence of gap to al¬ 
low a = 2A(0)/Tc > 3.53 (be. a values higher than 
the BCS value). The temperature dependence of each 
energy gap for this model can be approximated as®^: 
Ai(r) = Ai(0)tanh[1.82(1.018(^ - l))°-5i], where A(0) 
is the maximum gap value at T = 0. We adjust the tem¬ 
perature dependence of the London penetration depth by 
using the following expression: 

KbiO) ^Jo A(t,0) dE ^if2_A2(r,</.)’ 

(5) 

where A(T, <f) is the order parameter as functions of tem¬ 
perature and angle. For the two-gap model, is cal¬ 
culated as®”^: 

KbiT) = r\^\T) + (1 - r)X^\T), ( 6 ) 

where 0 < r < 1. Equations (5) and (6) are used to 
introduce the two gaps and their appropriate weights. 

The best description of the experimental data is ob¬ 
tained using values of Ai/ksTc = 1.72±0.3, 1.9±0.3 
and 1.74±0.25, Az/fcsTc = 0.49±0.3, 0.68±0.3 and 
0.59±0.25, and r = 0.2±0.1, 0.32±0.2, and 0.48±0.2 for x 
= 0.092, 0.096, and 0.1, respectively. The calculated pen¬ 
etration depth data are represented by the solid red lines 
in Fig. 5. It is noteworthy that our extracted gap values 
fit to the two-band s-wave fit for Ba(Fei_a;COa;) 2 As 2 ^^. 
Our investigated gap values for Hd and specific heat 
measurements have been found to be similar to those 
of values reported from the in-plane thermal conductiv- 
ity^®. In addition, our results are consistent with the 
BaFe 2 -a;Co 2 ;As 2 system, in which the superconducting 
energy gap does not contain a line of nodes anywhere on 
the Fermi surface, at any doping®®. On the other hand, 
the value of the gap amplitudes obtained for these SC 
samples scales relatively well with its Tc in light of the 
recent results for the Fe-based superconductors^®. Inter¬ 
estingly, one can notice that the extracted ratio for the 
anisotropic s-wave order parameter a is smaller than the 
BCS value, which points to the existence of the large gap. 


For the sake of comparison, we have summarized 
the Tc, the SDW transition T^r, 7 r, 7 „, the univer¬ 
sal parameter ACei/ynTj., and the values for the gaps 
Al, As for BaFei.908Nio.o92As2, BaFei,904Nio.o96As2, 
and BaFei. 9 oNio,ioAs 2 extracted from specific-heat and 
lower critical field (i7ci) measurements along with other 
hole-doped 122 materials in Table I. The A^/Ag 
ratio of the investigated systems in this work is 
found to be lower than in Bao, 6 Ko. 4 Fe 2 As 2 ^^ and 
Cao. 32 Nao. 68 Fe 2 As 2 ^^ systems, but this ratio is higher 
than the Bao, 65 Nao. 35 Fe 2 As 2 sample extracted from ear¬ 
lier specific heat measurements (Table I). The gap magni¬ 
tudes are scattered for different systems within the doped 
BaFe 2 As 2 . As mentioned above, the presence of a finite 
7 i. term is common in both electron- and hole-doped 122 
compounds. Most remarkably, assuming a SC volume 
fraction in our investigated SC samples (jn — lr)Hn ~ 
87.3, 96.1, 93.4% for x = 0.092, 0.096, and 0.10, respec¬ 
tively, which is in fair agreement with our magnetization 
data. Additionally, the relative weight of each contribu¬ 
tions illustrates that 72 / 7 ™ is always larger than 71 / 7 n 
indicating that the major gap develops around the Fermi 
surface sheet that exhibits the largest DOS. Theoreti¬ 
cally, in a two-band model that oc is expected 

in the interband weak-coupling limit®®. 

It is interesting to compare the present results with the 
other works for the most studied 122-based superconduc¬ 
tors in which the electron pairing mechanisms are still 
fairly under debate. For hole-doped Bai_a;Ka;Fe 2 As 2 , 
heat-transport measurements have claimed the possi¬ 
bility of line nodes in the SC gap in the underdoped 
regime®^. The similar nodal gap has also been ob¬ 
served in the heavily hole-overdoped Bao.iKo. 9 Fe 2 As 2 
and KFe 2 As 2 ^®’^^’®^. Interestingly, the isovalent substi¬ 
tution in Ba(Feo,64R.uo.36)2As2 and BaFe2(Aso,67Po.33)2 
showed a large residual in thermal conductivity and \/1T 
dependence, evidencing the presence of nodes in the SC 
gap®®. For electron-doped systems similar to the cur¬ 
rent study, the field dependence of the specific heat of 
both underdoped and overdoped Ba(Fei_a;Co 2 ,) 2 As 2 ex¬ 
hibits a Volovik-like nonlinear behavior, indicative of 
nodes in the SC gap®^, while nodeless gaps have been re¬ 
ported in underdoped compounds®®. Penetration depth 
experiments with a careful analysis of the SC state on 
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Ba(Fei_xCoa;) 2 As 2 concluded the possibility of node¬ 
less and nodes in the SC gap depending on the dop¬ 
ing level®®’®®. Very recently, the SC gap structure of 
(Bai_a;Ka;)Fe 2 As 2 was observed to vary with the compo¬ 
sition from two nodeless isotropic SC gaps at the optimal 
doping to a strongly anisotropic gaps at the end of the 
SC dome at a; = 0.16 [37,67]. In addition, the superfluid 
density of Ki_ 2 ,Naa;Fe 2 As 2 in the full temperature range 
follows a simple clean and dirty d-wave dependence, for 
pure and substituted samples, respectively®®. Near opti¬ 
mal doping for both hole- and electron-doped 122 com¬ 
pounds, various experiments have clearly demonstrated 
multiple nodeless SC gaps®’^^’^®’^^. In fact, it is hard to 
get a simple pairing mechanism from such complex situa¬ 
tion of the SC gap. The two independent techniques used 
here provide the self-consistent and convincing evidence 
for the nodeless gap in BaFe 2 -a;Nia;As 2 superconductors 
covering the underdoped to the optimal doping. 

Although in the current work we presented self- 
consistent data obtained from both magnetic penetration 
depth and specific heat measurements, some theoretical 
and other experiments also suggest a complicated pair 
symmetry for most iron-based superconductors, includ¬ 
ing various scenarios as mentioned above. However, it 
is important to emphasize that our investigated systems 
near optimal doping definitely underly and are consistent 
with nodeless multi-gaps in iron-arsenide multiband su¬ 
perconductivity in the presence of SDW, probably in the 
weak coupling regime. 

IV. CONCLUSION 

To summarize, from an extensive thermodynamic 
study of high-quality BaFe 2 -a;Nia;As 2 single crystals we 
have found that the magnetization loops exhibit a second 


peak, which is pronounced up to temperatures close to 
Tc- The main results are as follows, (i) Using the specific 
heat of a non-SC sample BaFei .75Nio. 25 AS 2 as a reference, 
we are able to separate the electronic specific heat from 
the phonon contribution for the SC samples down to T = 
0.4 K. (ii) Both the normal-state Sommerfeld coefficient 
and the jump of the specific heat AC/Tc are found to 
increase with Ni doping, indicating the strong competi¬ 
tion between superconductivity and magnetism, (iii) For 
all our SC samples, the electronic specific heat displays 
a pronounced anomaly at and a small residual part 
at low temperatures in the SC state, (iv) The observed 
temperature dependencies of C^i/T and are incon¬ 
sistent with a single BCS gap as well as with a d-wave 
symmetry of the SC energy gap. Instead, our analysis is 
consistent with the presence of two s-wave-like gaps in 
the nearly optimally doped compounds. 
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